Glassy properties of the Kawasaki dynamics of two-dimensional ferromagnets 
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We study numerically the Kawasaki dynamics of the 2d Ising model. At large time we recover 
the coarsening growth as l c (t) oc t 1 ^ 3 . At shorter time however, the system enters a metastable 
glassy regime that displays an extremely slow growth and non-trivial violations of the fluctuation 
dissipation theorem similar to those observed in spin glasses: this is one of the simplest system where 
such violations occur. We also consider Potts models, where a similar behavior is observed, and 
the model of Shore and Sethna where the domain growth is also slow, but where violations of the 
fluctuation dissipation theorem are trivial. We finally comment on these violations in the context of 
activated coarsening, and on similarities and differences with the glass transition phenomenology. 
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Last decades have witnessed many progresses in the 
study of out-of-equilibrium dynamics and coarsening, or 
domain growth, processes Q. In many simple models, 
the dynamics after a quench to the ordered phase is un- 
derstood in terms of domains growing with time as t 1 / 2 
with a Glauber (or non conserved magnetization) dynam- 
ics, and as i 1 / 3 with Kawasaki Q (or conserved magne- 
tization) dynamics 0. The 2d Ising model plays an im- 
portant role in the field as one of the simplest models for 
coarsening and phase separation. Recently, motivated 
by spin glass theory, many studies tried to go beyond 
these considerations to investigate two-time quantities 
such as correlation and response functions, the relation 
between them, and the way the dynamic Fluctuation- 
Dissipation Theorem (FDT) is violated in the out-of- 
equilibrium regime |j- Lll LJ ■ A generalization of the FDT 
can be formulated and amounts, within the mean field 
theory of aging, to the introduction of an effective tem- 
perature T e ff, it is widely accepted that T e ff — oo for 
coarsening process jfj , as confirmed by simulations in the 
Ising model with Glauber dynamics Q. In glasses and 
spin glasses however, violations of the FDT are highly 
non trivial |3|,y,|8|; finding lattice models with such prop- 
erties has recently triggered a lot of attention 0, [lfj ■ 

In this letter, we study the 2c? Kawasaki Ising model_ 
Even if the domains grow as i 1 / 3 at large time [ill 
we feel interesting to study the short time behavior as 
pre- as ymp totic regimes are already non trivial in the Id 
chain |l3|, We will demonstrate the existence of a 
glassy early regime, where domain growth is extremely 
slow, that displays a violation of the FDT similar to those 
of more complex frustrated models such as spin glasses 
or glass models. At large time we recover the Glauber- 
like behavior. The same phenomena arises with Potts 
variables but not in the model of Shore and Sethna > 
although it displays a logarithmic growth. We finally 
discuss our results in the context of glassy dynamics. 

Model and methods - We consider the Ising model 



whose critical temperature is T c w 2.27. To increase 
the efficiency of the Monte-Carlo (MC) simulations, we 
use a few tricks: first, we use the method described 
in [lij ]. a generalization to Kawasaki dynamics of an al- 
gorithm 0, 0] allowing a computation of the linear re- 
sponse to a field without physically adding a magnetic 
field in the simulation. Many two-time quantities for dif- 
ferent waiting times were thus computed in the same run. 
Secondly, we speed up simulations using the continuous 
time method of instead of choosing at random a 

pair of neighboring spins and to try to exchange them at 
each MC step, we keep track of all the Nf, broken links in 
memory, choose one at random, increase the clock by the 
time needed for the system to find it (a Poissonian vari- 
able with mean 2L 2 /Ni,; time units being in MC steps 
per site) and finally update the value of the correspond- 
ing pair of spins using the Heat Bath method. Since both 
tricks can be used at the same time, our approach allows 
a very good numerical determination of correlation and 
response functions at long time, even for large (L « 10 3 ) 
systems. Finally, we determined the coarsening length 
l c (t) by computing the excess energy 0,0,0] with re- 
spect to equilibrium l c (t) = —E eq / (Eft) — E eq ) (where 



E eq is given by the exact solution 
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nition l c (t) is 0(1) in the high temperature phase) and 
checked that the length defined as the first zero of the 
spacial correlation function jjjj gives similar results. 

Growing length beyond the i 1 / 3 law - Let us first 
consider in details the domain growth process. It is now 
clear that it behaves as i 1 / 3 at late time but first nu- 
merical studies failed to find this exponent and, after ini- 
tial claims in favor of a logarithmic coarsening 19], Huse 
stressed the importance of activated phenomena and en- 
ergy barriers [llj: due to the spin exchange dynamics, 
some moves need activation before happening, which cre- 
ates different time scales. For instance, zero tempera- 
ture simulations converge toward highly non trivial con- 
figurations far from equilibrium jlflj ]. where they remain 
blocked forever and therefore, since first excitations cost 
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FIG. 1: Coarsening length versus time at low T in a log-linear 
plot. After a finite time (O(10 2 )) the dynamics reaches the 
zero temperature plateau. After an time T\ oc e 4/3 , it enters a 
regime with a slow domain growth and beyond time T2 oc e 8/3 
eventually enters the t 1 ' 3 Kawasaki scaling (n and T2 are 
indicated for T = 0.25 on the picture). In inset: effective 
growth exponent versus ? c (£) _1 , that converges to 1/3. 



an activation energy SE = 4, any finite temperature sim- 
ulations for a time shorter than t\ = will behave as 
those at T = 0. Going at time larger than rj is not even 
enough to access the late time regime and in fact it is 
only when t ^> t-i = e s @ [ill Il2| that the system reaches 
the proper power-law behavior. 

Doing long quenches at many low temperatures, we 
investigate precisely the behavior of l c (t) and summarized 
our results in FIGQwhere three different regimes for the 
dynamics after a quench to T <C T c can be observed: 

(A) the Zero Temperature regime: after a finite time, one 
reaches the plateau characteristic of the T — dynamics, 

(B) the Glassy regime: after an activation time t\ oc e 4/3 , 
the dynamics leaves the plateau and enters a regime with 
a very slow domain growth. Actually, l c (t) seems even 
su b- log arithmic at large time, suggesting, as first intuited 
by [llj, that l c (t) saturates and reaches a second plateau 
at very large t (this is however quite invisible at the time 
scales considered here), and (C) the Asymptotic regime: 
after a time T2 oc e s/3 , the dynamics enters the canonical 
regime where l c (t) ~ A + Bt 1 ^ 3 , as can be seen from the 
behavior of its effective growth exponent (inset in FIG^I. 

This glassy regime, whose name will be justified later, 
thus corresponds to a metastable state with lifetime 
0(t2) that lasts extremely long at low temperature. 
That l c (t) seems to saturate indicates that the system 
is growing small organized domains that then remain 
blocked because further growth would need a move cost- 
ing SE — 8 (i.e. flipping a spin on domain edges). This 
looks very similar to the polycrystalline picture (a collec- 
tion of small crystalline zone) observed in the early stages 
of the dynamics of some model glass . We also noticed 
that the way the system escapes from these metastable 
states may involve cooperative phenomena and non just 
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FIG. 2: Up (2.a) : FDR plot at T = 0.7T C . At long time, 
curves converge slowly to the expected limit, as for Glauber 
dynamics However, at short time, curves are shifting up as 
t w increases, like in spin glass models. Down (2.b): FDT for 
T — 0.25T C , in the glassy regime only; curves are superposing 
very well on two decades in t w , and correspond to a non trivial 
quasi-equilibrium regime. In inset the corresponding X(G). 



a single spin flip of energy 8: modifying the dynamics in 
such a way that only the moves that cost 8E < 4 are 
allowed, we found that the dynamics still converge to the 
t 1 / 3 regime after a time, although a bit larger, still scaling 
as e s @ . To emphasize the differences between the dynam- 
ics in this glassy regime and the late time coarsening, we 
turn now to FDT violations. 

Violation of the FDT — Studies of the Fluctuation- 
Dissipation Ratio (FDR) 0, Q are based on the com- 
parison of how spontaneous and induced fluctuations do 
relax: one measures the two-time correlation function 
C(t,s) — (S(t)S(s)} and the associated response to a field 
function R{t, s) = d{S(t))/dh(s), and defines the FDR 
X(t, s) through the formula T R(t, s) = X(t, s) d s C{t, s). 
At equilibrium the FDT holds, thus X = 1. In the off- 
equilibrium regime however, the FDT is violated. In the 
large time limit (s,t — ► oo with C(t,s) — > q) the FDR 
X(t, s) converges to a limiting function X(q) whose phys- 
ical meaning has been derived within the mean field the- 
ory of aging, where it has been shown 5] that it can 
be computed using the overlap pdf P(q) in the thresh- 
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old states (the states reached by the dynamics on very 
large time 0) using X(q) = x(q) = f£ P(q')dq'. The 
interpretation of X(q) in finite d models, the meaning 
or even the existence of a unique effective temperature 
T e ff — T/X play a central role in present research in 
off-equilibrium physics U Experiments and simulations 
consider integrated quantities, in which case one gets 



x(t,t w ) = 



1-C(t,t w ) 



T, 



eff 



(1) 



where \(t,t w ) is the linear susceptibility at time t that 
results from a magnetic field switched on from times t w 
to t. A popular presentation of such data is a plot of 
Tx(t, t w ) versus the C(t 7 t w ): one finds that FDT holds at 
short time, when C(t,t w ) is large (in the so-called quasi- 
equilibrium regime) so that T e ff = T until C(t, t w ) drops 
and reaches a given value qEA Lil (equal to the exactly 
known equilibrium magnetization in the Ising model with 
Glauber dynamics) and FDT is violated for C(t,t w ) < 
qEA where the introduction of T e ff > T is thus needed. 

It is widely believed that any coarsening processes is 
characterized by T e tf — oo. This is quite expected since 
the P(q) corresponding to ferromagnetic equilibrium is 
trivial. Here, however, there is the early regime where 
the dynamics is reaching metastable states far from the 
equilibrium ones, so, inspired by the mean field interpre- 
tation of FDT violations, we expect that in the glassy 
regime they are not ruled by the equilibrium ferromag- 
netic P(q) but by the glassy metastable states so that the 
FDR is non trivial. We present our data for T = Q.7T C in 
FIG|21a: at large time, the curves converge very slowly 
toward the expected limit, approaching it downward as 
in the Glauber case (a direct comparison with 0, H| is in- 
structive). This excess response signal can be understood 
by noticing that spins on domain walls have low local 
fields and thus participate a lot to the global magnetic 
response. At short time l c (t) is quite low so there is a 
large excess density of domains that eventually decreases 
at larger time as domains annihilate: this gives rise to the 
bump in the FDT plot and explain the slow downward 
convergence of the susceptibility as t w increases. All this 
is similar with what happens with Glauber dynamics • 

However, a new phenomenon arises at shorter time: in 
FIG|21a curves are first drifting up as the waiting time 
t w increases. This is in sharp contrast with what is ob- 
served with Glauber dynamics and quite similar to spin 
glass behavior 8]: in this regime the system is some- 
how creating more response as t w increases. These phe- 
nomena arise for time scale 0(r 2 ) « I0 2 ; this suggests 
that, staying long enough in the glassy regime, by work- 
ing at a lower temperature, we should observe clearer 
FDT violations. This is indeed the case (see FIG|21b) at 
T = 0.25T C (where t-i rs 10 6 ) where we observe a super- 
position of FDR curves on more than two decades and 
obtain a plot very similar to those of complex systems 
as glasses and spin glasses. This shows that, between 
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FIG. 3: FDR plot for the 3d Shore-Sethna model with 
competitive ferromagnetic (J = 1) I s * neighbor and anti- 
ferromagnetic (J — —1/6) 2 nd neighbor interactions. We ob- 
served a very slow domain growth (inset) but a trivial FDR. 



times T\ and T2, the dynamics of the Kawasaki model 
in the glassy regime is very sensitive to the metastable 
states and, as expected, exhibits a non trivial FDR. It 
is interesting that ideas arising in mean field disordered 
complex systems are realized in a certain time range in 
the 2d Ising model, which as opposed to other simple 
models where FDT violations where observed [3, lift |2fJ , 
is not critical or close to a critical point, not disordered 
and not even statically frustrated. Despite that, it has 
the striking properties of behaving like a ferromagnet at 
long time but as a spin glass at short time, suggesting 
that any coarsening model with similar activated regimes 
could behave like a glass at short times. 

Logarithmic coarsening — To investigate the general- 
ity of our results, we followed Shore and Sethna fSSlflif 
who introduced a 3d Ising model with ferromagnetic first 
neighbors but weak anti-ferromagnetic second neighbors 
interactions, and argued that it has barriers growing with 
the size l c (t) of the ferromagnetic domains, thus leading 
to an activated logarithmic coarsening . We repeated 
our FDT study in the range of temperature and param- 
eters where SS claimed to observe a logarithmic growth 
and computed l c (t) using the same definition as [l5|- We 
present our results in FIGE3 although we find data com- 
patible with a logarithmic l c {t) at large time, x 1S very 
small, the FDR seems trivial and all our FDT plots look 
like what is obtained for usual coarsening process 0; 
in particular all curves are shifting downward as t w in- 
creases, as opposed to FIGEJa. This demonstrates that it 
is not just the slowly growing length that creates the non 
trivial FDR in FIG[21 There, the dynamics clearly dif- 
fers from usual coarsening, even logarithmic, and reaches 
some non trivial quasi-cquilibrium states that seem sta- 
ble until times of 0(t2); this is in contrast with the log- 
arithmic coarsening of the SS model which is just going 
slowly toward the ferromagnetic state, and where no such 
equilibration in a metastable state could be defined. 
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FIG. 4: FD in the g-Potts model for q = 10 (inset, T = 0.8T C ) 
and q = 15 at T — 0.5T C . As for Ising, one sees a crossover 
from the early regime to coarsening for t » T2 ~ O(10). 



Potts variables and glasses - The present study is 
easily extended to the q-states Potts model [2]) (where 
T c = q/ log (1 + y/q)) and we checked that a similar phe- 
nomenology for the FDR and the slow dynamics arises: 
the same activated processes give rise to glassy behavior 
(see FIG0J. This is interesting as Potts models have a 
I s * order transition for q > 4, hence we obtain a model 
with many features of the glass phenomenology such 
as FDT violations, melting transition and metastable 
states. Moreover, like other models with much more com- 
plicated Hamiltonians recently developed to study ana- 
lytically the glass transition 22] , our Ising or Potts mag- 
nets with Kawasaki dynamics have a non trivial mean- 
field limit on Bethe lattices, where they display a spin 
glass behavior (within the cavity method [23], the spin 
glass solution applies to ferromagnetic couplings when 
the magnetization is fixed to zero) so we may ask if this 
model is pertinent as a model glass and how it compares 
with ? As in any realistic models, the low T phase 
eventually crystalizes at large time once the free energy 
barriers to crystal nucleation are crossed (if the nucle- 
ation time is shorter than the equilibration time of the 
metastable glass state, the later cannot be thermodynam- 
ically defined: this is the solution Kautzmann originally 
proposed to its own paradox (2^|). A lifetime of 0(e 8/3 ) 
being not very large close to T c , it is quite easy in this 
model to overcome these barriers in a cooling experiment 
so that no supercooled state and no glass transition can 
be observed (this is also what happens in the model of 
obtained with a very frustrated plaquette Hamiltonian) : 
here, it is only by quenching the system fast enough to 
low T that it shows its glassy nature. Precautions are 
thus needed when using mean field methods; it suggests 
however that a similar model with barriers large enough 
not be overcome in the time scale of the simulation would 
be a good model; this is probably what happens in mod- 
els [22| which are believed to be good glass formers. 

Discussion — We studied the off-equilibrium dynam- 
ics of 2d Kawasaki magnets and demonstrated the ex- 



istence of a early extremly slow glassy regime with a 
non trivial violation of the FDT similar to those of spin 
glasses. We observed the same phenomenology with 
Potts variables but not in the SS model, pointing out 
the difference between this dynamical behavior and the 
usual coarsening, however slow it may be. We believe 
that this behavior is quite generic and should be observed 
(as it has already been in some specific cases 0,0]) in 
the early stage of other coarsening model (provided that 
they share similar activation properties) or in some lat- 
tice model for glasses. Since non trivial FDT violations 
have been observed in many systems it is interesting 
that a simple model such as the 2d Ising ferromagnet, 
whose statics solution is analytically known for decades, 
has a dynamics that already shares many properties with 
those of much more complex materials: very few ingre- 
dients are needed for that. Finally our results show, as 
was pointed out in [2(J , that great cares should be taken 
when using dynamical simulations to probe the statics. 
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